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A new wavefunction which improves the Gutzwiller-type local ansatz method has been pro- 
posed to describe the correlated electron system. The ground-state energy, double occupation 
number, momentum distribution function, and quasiparticle weight have been calculated for 
the half-filled band Hubbard model in infinite dimensions. It is shown that the new wavefunc- 
tion improves the local-ansatz approach (LA) proposed by StoUhoff and Fulde. Especially, 
calculated momentum distribution functions show a reasonable momentum dependence. The 
result qualitatively differs from those obtained by the LA and the Gutzwiller wavefunction. 
Furthermore, the present approach combined with the projection operator method CPA is 
shown to describe quantitatively the excitation spectra in the insulator regime as well as the 
critical Coulomb interactions for a gap formation in infinite dimensions. 

KEYWORDS: variational method, electron correlations, Gutzwiller wavefunction, local ansatz, 
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1. Introduction 

Variational method has been a useful tool to investigate the ground-state properties of 
correlated electrons from molecules to solids over half a century. In the method, a minimum 
basis set to describe correlated electrons is constructed by applying one-particle, two-particle, 
and higher-order particle operators onto the Hartree-Fock wavefunction, and their amplitudes 
are chosen to be best on the basis of the variational principle. Among various trial wavefunc- 
tions, the Gutzwiller wavefunction is one of the simplest and popular wavefunctions in solids. 
It was introduced by Gutzwiller to clarify the role of electron correlations in metallic ferromag- 
netism.^'^ The idea is to reduce the amplitudes of doubly occupied states on the local orbitals 
in the Hartree-Fock wavefunction by making use of a projection operator nj(l — gnii^ni[). 
Here riio- is the number operator for electrons on site i with spin a. The variational parameter 
g is determined by the minimization of the ground-state energy. Later, it has been pointed 
out by Brinkman and Rice that the Gutzwiller wavefunction describes the metal-insulator 
transition.^ Because the Gutzwiller method is a nonperturbative approach, it has extensively 
been applied to the strongly correlated electron systems such as the heavy fermions, high-Tc 
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cupratcs, and other transition metal oxides.^ 

Although the Gutzwiller ansatz approach (GA) captures the physics of electron correla- 
tions and is useful for correlation problems, it was not so easy to apply the method to realistic 
Hamiltonians. Stollhoff and Fulde proposed an alternative method called the local-ansatz 
approach (LA), which is simpler in treatment and applicable to realistic Hamiltonians.^'^ 
The LA takes into account the excited states created by local two-particle operators such as 
{Oi} = {(5nj|(5njj}, and determines their amplitudes variationally. Here Srii^j = rii^ — {nia-)^, 
{'nia)o being the average electron number on site i with spin a in the Hartree-Fock approxima- 
tion. The LA has been applied to many systems such as molecules, polyacetylene, transition 
metals, semiconductors, and transition metal oxides.''' 

The LA is useful for understanding correlation effects in actual materials. The application 
however has been limited to the weakly correlated region because of the difficulty in evaluation 
of the higher-order terms in average quantities. In the weak interaction limit, the Hilbert 
space expanded by the local operators is, however, not enough to describe exactly the weakly 
correlated region; the LA does not reduce to the second-order perturbation theory in the weak 
correlation limit. The same difficulty also arises in the original Gutzwiller wavefunction even in 
infinite dimensions. In the present paper, we aim to solve these problems in the LA introducing 
a new wavefunction with momentum-dependent variational parameters, and demonstrate that 
the new approach much improves the LA in the weak and intermediate correlation regimes. 
In the following, we call the new approach the MLA (the LA with momentum dependent 
variational parameters). 

We write down our wavefunction for the single-band Hubbard model in §2. The idea is to 
choose the best local basis set obtained from the two-particle excited states in the momen- 
tum representation by projecting out those states onto the local subspace and by controlling 
the amplitudes of the excited states in the momentum space. We calculate the ground-state 
energy within a single-site approximation. Using the variational principle, we determine the 
momentum-dependent variational parameters. The ground-state energy obtained from our 
wavefunction agrees with the result of the second-order perturbation theory in infinite dimen- 
sions in the weak interaction limit, and reduces to the correct atomic limit in the case of the 
half-filling. 

In §3, we present the results of numerical calculations for the half-filled band Hubbard 
model in infinite dimensions to examine the validity of the new wavefunction. We calculate the 
correlation energy, the double occupation number, the momentum distribution function, and 
the quasiparticle weight as a function of the Coulomb interaction energy parameter. We verify 
that the present approach improves both the LA and the GA in the weak and intermediate 
Coulomb interaction regimes. In particular, we demonstrate that the momentum distribution 
calculated by our wavefunction (the MLA) shows a distinct momentum dependence, and 
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is qualitatively different from those obtained by the LA and GA leading to the constant 
values of the distribution function below and above the Fermi level. In §4, we present an 
example of applications to excitation problems. We have recently developed a self-consistent 
method^ to calculate the excitation spectra from the retarded Green function by making use 
of the projection operator technique and the effective medium within the coherent potential 
approximation (CPA).^'^'^ The method called the projection operator method CPA (PM- 
CPA)^ is equivalent to the many-body CPA/"^ the dynamical CPA/"^' and the dynamical 
mean field theory /^"^^ and treats the dynamics and the static correlations separately in the 
calculations. We calculate here the excitation spectra combining our variational method with 
the PM-CPA. We show that the calculated spectra in the insulator regime quantitatively 
agree with the results of the numerical renormalization group (NRG) calculations.^" In the 
last section, we summalize our results and discuss future problems. 



2. Local Approach with Momentum Dependent Variational Parameters 

We consider in the present paper the single-band Hubbard model defined by 

icr ija i 

Here eo {h) is the atomic level (magnetic field), Uj is the transfer integral between sites i and j. 
U is the intra-atomic Coulomb energy parameter, aj^ (oia) denotes the creation (annihilation) 
operator for an electron on site i with spin a, and riia = a\^o,i(j is the electron density operator 
on site i for spin a. 

In the Hartree-Fock approximation, we replace the many-body Hamiltonian (1) with an 
effective Hamiltonian for independent electrons, 

Ho = ^ tija a\^aj^ - ^ {n^)o{nii)o , (2) 

ija i 

and approximate the ground-state wavefunction |*) with that of the Hartree-Fock Hamil- 
tonian i^o, «-e., I^o). Here tija = (eo + U{ni-a)o — crh)dij + tij{l — Sij). (~)o denotes the 
Hartree-Fock average (^z!>o|(~)|0o)) and {nia)o is the average electron number on site i with 
spin a. The Hamiltonian (1) is then expressed by a sum of the Hartree-Fock Hamiltonian and 
the residual interactions as 

H = HQ + uY,Oi . (3) 

i 

Here Oj = 5ni^5nii and 5ni^ = nia - {nia)o. 

In the local-ansatz approach (LA),^ we take into account the Hilbert space created by 
operation of the residual interaction {Oj} onto the Hartree-Fock state l^o). Introducing a 
variational parameter t^l^ into the basis set {Oj}, the LA wavefunction for the ground state 
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is written as 

|*la)= [n(l-^LAO«)]l'/'o) . (4) 

i 

The LA is different from the Gutzwiller ansatz wavcfunction I'I'ga) = Wii^ ~ Q'^^i^'^^il) |?^o) 
in which the doubly occupied states are exphcitly controlled by a variational parameter g, 
and simplify the evaluation of the physical quantities in the weakly correlated region. 

As we have emphasized in the introduction, the LA does not lead to the exact result in the 
small U limit because it makes use of a limited local subspace. In fact, the LA wavefunction 
(4) may be expanded formally in the weak interaction limit as 

|*la) = |<^o) + |'^i)la + --- , (5) 
|0i)la = -X1 ^ (^iN)(^l^i)(4K)(^l^2)?7LA^(alu"'=2i)^^4iT"'=iT)l^o) • 

i kik^k'zk'^ 



Here {i\k) = exp(— ifc • R.i)/\/N is an overlap integral between the localized orbital and the 
Bloch state with momentum fc, Ri denotes the atomic position, and N is the number of sites. 

(aka) denotes the creation (annihilation) operator for an electron with momentum k and 
spin a, and S{al,^aka) = al/^Ofca - (aI/^afea)o- 

The Rayleigh-Schrodinger perturbation theory, on the other hand, yields the following 
form 

l*) = l<^o) + |^i) + -- - , (7) 

i kik[k2k2 

Here /(e) is the Fermi distribution function at zero temperature, and = ^ka — At- is 
the Fermi level, eka is the Hartree-Fock one-electron energy eigen value given by eka = eo + 
U{ni-a)o + Cfe — (^h, and is the Fourier transform of Uj. 

Equation (8) compared with eq. (6) manifests that one has to take into account the 
momentum dependence of the variational parameters to improve the LA. We propose in the 
present paper the following wavefunction with momentum-dependent variational parameters 

{Vk'2k2k[ki}- 

\^) = Y[{i-di)\^o) , (10) 

i 

Oi= Yl (^iK)(^l^i)(^2K)(^l^2)%^fc2fcjfei5(aJ.,^afc2^)5(aJ.,^afc,t) . (11) 

kik2k[ k'2 
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The operator Oj is still localized on site i because of the projection {k[\i) {i\ki) {k2\i) {i\k2) ■ It 
should be noted that Oj 7^ Oj and OiOj 7^ OjOi {i 7^ j) in general. These properties do not 
cause any problem when we make a single-site approximation. When we treat the nonlocal 
correlations we have to adopt the symmetrized wavefunction in general. The wavefunction 1^*) 
reduces to I^'la) when {r]y_^k^k[kx\ become momentum-independent. 

The variational parameters are determined by minimizing the ground-state correlation 
energy Ec- 

= {H) - {H)o = . (12) 

Here H = H - {H)o. 

Calculation of the correlation energy with use of the new wavefunction is not easy in 
general. But, one can evaluate it within the single-site approximation. As shown in Appendix 
A, the average {A) of an operator A = A — {A)q with respect to the wavefunction (10) is given 
in the single-site approximation as 

a)=v«i^^Mi^. (13) 

i ((i-o;)(i-o,))o 

By making use of the above formula, one can obtain the correlation energy per atom. 

^ ^ -{OlH)o-{m)o + {OlH6i)o .^4. 
l + (OjO.)o 

Here we assumed that all the sites are equivalent to each other and we made use of the fact 

(oi)o = mo = 0. 

Each term in the correlation energy (14) can be calculated by making use of Wick's 
theorem as follows. 

kik2k[k'2 j 

x%iit2feifei/(efeiT)(l - /(4iT))/fei)(l - /(ejfc^)) , (15) 
{d}H)o = {Hd,)l , (16) 

(dj#0,)o = (^Wi){ki\i)ii\k'2){k2\i)ri\k2k',k, 



kik2k[k'2 

x/(6fc,T)(l - /(4;T))/(e~fe2i)(i - E {k3\i){i\h){k',\i){i\k,) 

/i<3 ^4 ^3 ^4 

X (^^Ek'2k2k[ki^kik3^k{k'Jk2k4^k'2k'4^ + U^'^k^k'^^kik'^kik'^ka^ Vk'Juk'^ks i (17) 

Uk'2k2k[kik'^k4k'sk3 = U Y[{j\ki){k3\j)fiek3t)h{k's - {k'i\j){j\k3){i - f(.h'3^))Skik3] 

j 

x[{j\k2){k4\j)f{ek4l)Sk'ki - (4|j)01fe4)(l - fiek'Am2kA , (18) 
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{0\0i)o = 4 E l%ifeikifejV(efeiT)(l-/(4iT))/(efei)(l-/(4u)) • (19) 



Here AE'^.^^.^fc^y^j = e^/^ — e^jj, + e^'^i — efcjt is a two-particle excitation energy. 

The above expressions (15) and (18) contain nonlocal terms in the summation over j (i.e., 
"Ylij)- thus make additional single-site approximation called the = approximation.^^ 
In eq. (15), for example, we have 

Y.^mmi){mm2){ki\j)m){k2\j)m) = ^y.^^c^'^'"'"''''''^^'''''''^ ■ (20) 

The R = approximation only takes into account the local term (j = i) in the above sum- 
it /nt , 



mation, so that {HOi)o{= {OJH)q), and {OjHOi)o reduce as follows. 



(21) 



kik2k[k'2 



kikik^k'^ 
^-^fc^tefcjfci Vk'^kik^ki 



+ E /(^feT)/(^fc4i) Vk'2k4k[k3 - E /(^'=3T)(1 ~ /(efc^i)) Vk'^k2k[k3 



kski 



ksk' 



- ~ /(e^T))/(^fc4i)%^fc4fc^fci + ~ /(^^t))(1 - /(efc4i))%4fc2^fci} • (22) 

Variational parameters {r]k'2k2k[ki} ^'^^ obtained by minimizing the correlation energy e^, 
i.e., eq. (14) with eqs. (19), (21), and (22). The self-consistent equations for {^^fcjfcifci} 
single-site approximation are given as follows. 

{^^k!2k2k[ki ~ ec)%^fe2it'jifci 

+ [ E fi^kadfi^kiOVk'^kik'^ks - E /(^fc3T)(l ~ fihil))Vk'^k2k[k3 



kski 



kzk'. 



- E(l ~ /(efc^T))/(^fc4i)%^fe4^iki + ~ /(4^t))(1 - f{hii))r]k'^k2k'^ki = U . (23) 

Aig A/4: ^2 

Note that r]k'2k2k[ki should vanish when U ^ 0. Thus in the weak U limit, one can omit the 
second term at the l.h.s. (left-hand-side) of eq. (23). We then obtain the solution in the weak 
U limit as 



_ U 

^^k'2k2k[ki 



(24) 



In the atomic limit, AEj^/^f^^k^ki = 0- We find then a ^-independent solution being identical 
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2{OiH)o{Of)o 



(25) 



It is not easy to find the solution of eq. (23) for the intermediate strength of Coulomb 
interaction U. We therefore consider an approximate solution which interpolates between the 
weak and the atomic limits. Note that the second term at the l.h.s. of eq. (23) do not affect 
the solution in the weakly correlated limit as we have mentioned. Therefore we approximate 
{rjk^k^k'j^ki} the second term with the momentum- independent parameter rj which is suitable 
for the atomic region. We have then an approximate solution as follows. 

?7[l-r/(l-2(nit)o)(l-2(nii)o)] 



'^Ek'2k2k[ki 



(26) 



The best value of r] should be determined variationally, but wc make use of that in the LA for 
simplicity. Furthermore, we approximate ec in the denominator of eq. (26) with the correlation 
energy in the LA. 

Substituting the variational parameters (26) into eq. (14), we obtain the ground-state 
correlation energy. The each element in the energy is given as follows. 

{HOih = {0\h)1 

= C/2[l-r;(l-2(n,t)o)(l-2(n,;)o)] 

[Unden\p^iel)p^ie2)pi{e^)p^{e,)f{el){l - /(e2))/(e3)(l - /{e,)) 



I 



£4 — 63-1- 62 



{dlHdi)o = {olHodi)o + u{diOidi)o , 



(27) 
(28) 



{6lHo6i)o = [/2[l-r/(l-2(n,T)o)(l-2(n^)o)P 

Ylden P]{ei)p^ {62) Plies) Pi{e4) 



{olOidi)o 



x/(6i)(l-/(62))/(63)(l-/(64)) 



U^[l-r]{l-2{nii)o)il-2{nii)o)f 



( L - < + (2 - (1 
(64 — 63-1-62 — 61 — 60 



(29) 



Ylden Piiei)p^{€2)pi{e3)pi{e4 

n 

(ie5(ie6Pt(e5)Pi(e6)/(e5)/(e6) 



/(6i)(l-/(62))/(63)(l-/(64)) 



-/ 



64 — eg + 62 ~ €5 — ec 
(i65(i66Pt(e5)pj^(e6)(l - /(e5))/(e6) 

64 — 66 + 65 — 61 — 6c 



64 — 63 -I- 62 — 61 — 6c 

d65d66Pt(e5)p|(e6)/(e5)(l - fjee)) 
£6 - £3 + 62 - £5 - ec 

d65d66pt(65)p^(66)(l - /(e5))(l " /(cg)) 
^6 ~ ^3 -|- 65 — 61 — 6c 



,(30) 



f/2[l-r7(l-2(nit)o)(l-2(n^)o)]' 
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(31) 



(£4 - £3 + €2 - ei - ec)2 

Here Pa{e) is the density of states for the Hartree-Fock one-electron energy eigen values mea- 
sured from the Fermi level. 

Electron number (nj)(= X^(^(?^^^r)), the momentum distribution (nka), and the double 
occupation number (ni|ni|) are obtained from d{H)/de, d{H) /dika, and d{H)/dUi, respec- 
tively. Here ika = e^ — crh and is the Fourier transform of Uj. Making use of the single-site 
energy (14) and the Feynman-Hellmann theorem, we obtain the following expressions. 

{dinidi)Q 



{nil = {ni)o + 



{nka) = {nka)o + 



l + (Ola)o ' 



{OlOi)0 - {OiOi)o + {6lOiOi)o + Ea{ni-aMOjniaOi)o 



l + (OjOi)o 
Here Ui = Ui - {ni)o, hka = nka - {nk<7)o, and 

{dlhM)o = 2;72[l-77(l-2(n,|)o)(l-2(n^)o)]2 

P_,(ei)p_,(e2)pa(e3)p.(e4)p.(e5)/(ei)(l - /(e2))/(e3)(l - /(e4)) 



(32) 
(33) 
(34) 



X 



£4 — 63 + £2 — ei — ec 
1 - /(£5) /(£5) 



N{dlhkadi)o 



.£5 - £3 + £2 - €l - £c £4 - £5 + £2 - £1 - £c 

C/2[l-r/(l-2(n,T)o)(l-2(n,;)o)F 

dei(ie2(ie3/9_o-(£i)/9-a(£2)Po-(£3)/(£i)(l - /(£2))/(£3) 



(35) 



1 



fieka)) I ■ 



(ei - 62 + £fcCT - £3 - £c 



deide2(i£3P-o-(£i)/3-a(£2)Pa(£3)/(£i)(l - /(£2))(1 - /(£3)) 



(ei — £2 + £3 — ^ka — £c 



(36) 



{dlOi)o + {Oidi)o = 2C/[l-r/(l-2(n,t)o)(l-2(nii)o)] 

^ /- [nn^^n]pT(^l)PT(^2)Pi(£3)Pi(£4)/(£l)(l - /(£2))/(£3)(l - /(£4)) _ ^3^^ 
J £4 — 63 -I- 62 — £1 — £c 

It should be noted that the expressions of these physical quantities consist of the multiple 
integrals up to the 6-folds. One can reduce these integrals up to the 2-folds using the Laplace 
transform. Their expressions are given in Appendix B. 

3. Numerical Example: Half-Filled Hubbard Model 

We have performed the numerical calculations of the half-filled band Hubbard model in 
order to examine the properties of the local ansatz with momentum-dependent variational 
parameters (MLA). We consider here the hypercubic lattice in infinite dimensions, where the 
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1.8 




Fig. 1. Calculated elements {0\H)q (dashed curve), {0\HOi)Q (solid curve), and (0|0,)o (thin solid 
curve) as a function of the Coulomb interaction parameter U. Corresponding curves in the LA are 
shown by dotted curves. 



single-site approximation works best. The density of states for noninteracting system on the 
hypercubic lattice is given by p(e) = (1/y^) exp(— e^), in which the energy unit is chosen to 
be / dep{e)e^ = 1/2. 

Figure 1 shows the curves for {d\H)^, {dlHdi)o, and (0|a) in the correlation energy 
as a function of the Coulomb interaction energy parameter U. These are proportional to C/^ 
in the small U limit. The energy {0\h)q agrees with {0\HOi)Q in the small U limit as is 
verified from eqs. (27) and (28). This implies that ec ~ —{0\h)q in the small U limit. For 
large U , the energy {p\HOi)Q becomes smaller than {0\h)q because the former changes 
with {r]y_^k,^k[kx\ ™ ^ quadratic way, while the latter linearly depends on U as U r]y^i.^y^^^, 
and because {^^fcjfcifci} saturate when U is large. The renormalization contribution (0|0j)o 
tends to saturate with increasing U because of the same reason. As a consequence, the energy 
— {0\h)q forms the leading term in the correlation energy even for large U . 

The same quantities in the LA are obtained by the replacements T]i.i_^]^^y^j.^ — > '7lA' 
{d\H)o ^ VLA{OiH)o, {dlHdi)o ^ ril^{0,HO,)^, and (0|0,)o ^ r,lK{0^i)o in eqs. 
(27), (28), and (31). Their expressions are calculated analytically in the present case as 
r/LA(Oi^)o = ??LAf^/16, riip^iOiHOi)^ = vlJ^V^, viA{Of)o = r?iA/16, and 



1 11 
+ \ - + 



Vn \ IT 16 , . 

^LA = ij • (38) 

32 
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Fig. 2. The correlation energies vs. Coulomb interaction energy parameter U in the LA (dashed 
curve), the ML A (solid curve), and the GA (dotted curve). 



These quantities in the LA are also presented in Fig. 1 by dotted curves. The results in the 
LA describe those in the MLA rather well over a wide range of the Coulomb interaction. It 
should be noted that the curves in the LA deviate from those in the MLA even in the small 
U limit. A remarkable point is that (O|^)o is larger than that in the LA. This lowers the 
ground-state energy of the MLA. 

Figure 2 shows calculated correlation energy as a function of Coulomb interaction. The 
energy in the MLA is lower than that of the LA over all Coulomb interaction energy parameters 
U, verifying an improvement of the wavefunction. Moreover, the MLA wavefunction leads to 
the correlation energy lower than that of the original Gutzwiller Ansatz (GA) for the Coulomb 
interaction energy parameter U < U* = 3.28. We can expect that the MLA yields better 
results in the interaction range [0, U*], while the original GA should be better for U >U*. 

We present the double occupation number vs. Coulomb interaction curves in Fig. 3. The 
double occupancy in the uncorrected limit is 1/4, and decreases with increasing Coulomb 
interaction U. Both the LA and the MLA yield (nj-|-njj) = in the limit U = oo. The MLA 
suppresses {ni^riii) of the LA typically by about 10 % in the intermediate regime of Coulomb 
interaction. The double occupation number in the GA linearly decreases with increasing U 
and causes the metal-insulator transition at Uc2 = 8/-\/tt{= 4.51). The GA underestimates 
the double occupancy in the insulator regime because (rXj-fn^) should be finite even beyond 
Uc2 due to virtual exchange of electrons between the nearest neighbor atoms. Present result of 
the MLA indicates that the GA overestimates {ni^riii) for small U{< 2) and underestimates 
it at [/" ~ 3. 



10/23 



J. Phys. Soc. Jpn. 



Full Paper 



c 

«- 
c 



0.25 



0.2 



0.15 



0.1 



0.05 



LA 
MLA 
GA 



3 
U 



Fig. 3. The double occupation number vs. Coulomb interaction energy U curves in the LA (dashed 
curve), the MLA (solid curve), and the GA (dotted curve). 



The difference between the LA and the MLA is seen more clearly in the momentum- 
dependent quantities. Figure 4 shows the momentum distribution in various approximations. 
The distributions in the LA are constant below and above the Fermi level irrespective of U. 
This behavior is also found in the GA.^ The MLA shows a distinct momentum dependence 
of {nkc) via the energy e^. The results are in good agreement with those in the previous 
results of the RPT-1 (The first-order approximation in the renormalized perturbation theory) 
in the projection operator method CPA.^ The latter is exact up to the second order in U, and 
reproduce the Hubbard III approximation in the large U region. 

The jump at the Fermi level in the momentum distribution gives us the quasiparti- 
cle weight Z {i.e. the inverse effective mass). Calculated quasiparticle weight vs. Coulomb 
interaction curves are shown in Fig. 5. The quasiparticle weight in the LA changes as 
Z ={1- 3rilJl6)/{l + t/Ia/IS) and vanishes at C/c2(LA) = 24/V37r (= 7.82). In the GA,^ 
the quasiparticle weig ht changes as Z = 1 - {U/Uc2f. The curve in the GA agrees with the 
LA curve for small U. But it deviates from the LA when U becomes large, and vanishes at 
[/c2(GA) = 8/-v/7r(= 4.51). It should be noted that the GA curve strongly deviates from the 
curve in the NRG^'^ which is considered to be the best at present. The quasiparticle weight 
in the MLA much improves the LA; it is close to the curve in the NRG up to ?7 2.5, and 
vanishes at ?7c2(MLA) = 3.21. The latter should be compared with [/c2(NRG) = 4.10. We 
note that the wavefunction itself does not show the metal-insulator transition at Uc2 in the 
present approximation because the approximate expression of variational parameters (26) has 
no singularity at finite value of U. In this sense, the calculated Z and wavefunction are not 
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self-consistent in the present approximation. The values of Z obtained by the LA and MLA 
should be regarded as an estimate from the metallic side. 

4. Application to Excitation Spectra 

We demonstrate in this section that the present approach is also useful for understanding 
the correlation effects on excitation spectra. This can be made by combining the MLA with the 
projection operator method (see, for example, Chap. 6 in Ref. 7). The projection technique 
treats the dynamics of electrons and the static average separately, and the latter can be 
calculated by using the wavefunction method. 

We adopt again the half-filled band Hubbard model in infinite dimensions, and apply 
the projection operator CPA method (PM-CPA).^ In this methods, we describe the single- 
particle excitations by means of the Fourier transform of the retarded Green function, (aJ^^Kz— 
V)~^a}-^. Here the Liouville operator L defined by LA = [H, A]^ for an operator A describes 
the dynamics of electrons, z = uj + id, S is an infinitesimal positive number, and the inner 
product between the operators A and B is defined by {A\B) = {[A^B]^). We approximate 
in the PM-CPA the operator L by an energy dependent Liouvillean L{z) for an effective 
Hamiltonian with a coherent potential ^(z), ^i^)i^ia-'r^ija ^ij ^la^Jf^- "^^^ Green function 
for L{z) is given by 

H^) = / ^If^— . (39) 

J z — T>{z) — e 

Here p{e) is the density of states (DOS) per site for the noninteracting system. 
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Fig. 5. Quasiparticle-weight vs. Coulomb interaction curves in various theories. The RPT-1: dashed 
curve, the NRG: thin sohd curve, the LA: dotted curve, the MLA: solid curve, and the GA: 
dot-dashed curve. ^ 



To obtain the coherent potential ^{z), we consider an impurity system with Coulomb 
interaction on a site embedded in the coherent potential. The site-diagonal impurity Green 
function for the system is then given by 

G» (z) = - A» (z) + ~^ . (40) 

The self-energy M''\z) is expressed as follows for the half-filled band according to the RPT 
(Renormalized Perturbation Theory).^ 

A«(^) = ^?^^^, (41) 
1 + 4S(^)gJ^(z) 

G^\z) = {Al\{z - Lo{z) - Lf{z)Q)-'Al) . (42) 

Here A^^ = al^Srii-cr is an atomic operator expanded by the intra-atomic Coulomb interaction. 
Lq{z) = QL{z)Q and L^i\z) = QL^j\z)Q are respectively the coherent Liouville operator 
and the interaction Liouville operator in which Lf\z) is defined by L^p{z)A = [USni^Sriii — 

S(2;)njo-,^]- Operator Q (Q) denotes a projection operator which removes the original 
operator space {|ajj^)} (the atomic operator space {|^|(^)} )■ Note that the above expression 
(41) is exact, and Gq {z) denotes a screened memory function in which the dynamics for the 
strong atomic excitations have been removed. 

We consider here the lowest-order approximation (RPT-0)^ ; we neglect in Gq (z) a 'weak' 
interaction Liouvillean Lj {z)Q. Expanding the operator A'-^ by means of the nonlocal oper- 
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ators {al^6{al,_^ak"-o-)} as J2kk'k" (^ka^i^i'-a^k"-a){k\i){k'\i){i\k"), we reach a simple form 
of the screened memory function as follows (see eq. (74) in Ref. 8). 

^0) ^ f dede'de"p{e)p{e')pie")Xie,e',e") 
° ^ ^ J z- t{z) -e-e' + e" ' ^ ' 

Here we have omitted the spin dependence for simplicity. The correlation function X(e, e' , e") 
is given by 

X{ek,ek',ek") = e'(''-'"-'i+'")-'''(5(4,_^afe._,)5(4,,_^afe,_,)) 



kik'l 

The function X should depend on the momentum only via e^, an one-electron eigenvalue in 
the Hartree-Fock approximation. 

The coherent potential ^(z) is obtained from a self-consistent condition (i.e., the CPA 
equation^' ^^); 

G^^{z)=F{z). (45) 

Note that eq. (43) reduces to the second order self-energy when X{e,e',e") is treated by 
the Hartree-Fock approximation. Thus, the self-energy (41) yields the exact weak Coulomb 
interaction limit. The self-energy (41) also becomes exact in the atomic limit. 

In order to obtain the explicit expression for the self-energy from eq. (44), we adopted in 
our previous paper^ the Hartree-Fock wave function. We adopt here the new wavefunction 
(10) and the single-site approximation {i.e., R = approximation^^). We have calculated the 
function X(efe, e^', e^//) in eq. (43). Actual expressions used in the numerical calculations are 
given in Appendix C. The RPT-0 memory function obtained from the R = approximation in 
general does not satisfy the Fermi liquid condition. Therefore, the lowest order approximation 
is not applicable for the metallic state. Here we limit ourselves in the insulating state to 
demonstrate the quantitative aspect of the MLA within the lowest order approximation (RPT- 
0). 

Calculated excitation spectra in the insulating regime are presented in Fig. 6. The result of 
the MLA is compared with the Hartree-Fock one and the NRG which is considered to be the 
best at zero temperature. The spectrum with use of the Hartree-Fock wavefunction indicates 
the insulator. But the upper and lower band widths are broader than those of the NRG. Static 
correlations localize the electrons and suppress such band broadening. Resulting spectrum in 
the MLA reproduces well the NRG one. The agreement of the spectra implies a quantitative 
description of the site-diagonal Green function according to the Lehmann representation of the 
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Fig. 6. Single-particle excitation spectra calculated by means of the RPT-0 plus MLA (MLA-RPTO: 
solid curve), the RPT-0 plus Hartree-Fock approximation (HF-RPTO: dotted curve), ^ and the 
numerical renormalization group method (NRG: dashed curve). 



Green function G*^*) (z). In infinite dimensions, this means that the self-energy T.{z) is described 
quantitatively by means of the present theory because of the CPA equation G^''\z) = F{z) 
and eq. (39). 

The critical Coulomb interaction Ud for gap formation is obtained from a condition that 
the insulator solution lml](0+) = — oo disappears. This is equivalent to the following condition 
in the RPT-0 (see eq. (86) in Ref. 8). 

U = 4v^ . (46) 

Here C2 is the second moment of the memory function: C2 = / dede'de" p{e) p(e') p{e") (e + e' — 
e")^X{e,e',e"). 

For the wavefunction (10), the second moment C2 is given by the Hartree-Fock contribution 
02'^ = 3/8 + 3q;^/2 and the correlation correction as 

(2) 

C2 = + . (47) 

Here a = I/v/tt for the hypercubic lattice and a = 4\/2/37r for the Bethe lattice. The corre- 
lation contribution 03^^ in the second term at the r.h.s. of eq. (47) is given as follows. 

= 6 J (ieide2cie3p(ei)/9(e2)p(e3)(ei - €2 - e3)^[/(ei)/(e2)Ki(-ei - €2) 

-/(ei)(l - /(e2))AC2(e2, ei) + /(ei)(/(e3) - /(e2))A(-ei)] , (48) 
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Table I. Critical Coulomb interaction Ud for the Bcthc lattice in various approximations. MLA- 
RPTO: The MLA plus the lowest-order RPT in the PM-CPA (Present resuh), IPT: Iterative per- 
turbation theory,^^ ED: Exact diagonalization method, PSCT: Projective self-consistent tech- 
nique,^^ LMA: Local moment approach, NRG: Numerical renormalization group approach,^° 
1/U Exp.: 1/U expansion method.^^ 



MLA-RPTO 


IPT 


ED 


PSCT 


LMA 


NRG 


1/U Exp. 


3.36 


3.67 


3.04 


3.39 


3.41 


3.54 


2.97 



Ki(ejfc) = j dede'p{e)p{e')f{e){l - /(e'))Ko(eifc, e' - e) , (49) 

i^2iek, ^k') = J dedeV(e)p(e')/(e)(l - /(e'))Ko(eik - e, e' - e^/) , (50) 

A(efe) = J dep{e)il - /(e))Ao(efe, e) , (51) 

«o(efe, ek') = J dede'p{e)p{e')f{e){l - /(eO)r?(e' - e + ek)il{e' - e + e^O , (52) 

Ao(efe,efeO = j dede'de"p{e)p{e')p{e")f{e){l-f{e'))f{e") 

xr]{e' -e-e" + SkHe' - e - e" + e^O • (53) 

Here rj{ek) = U/{ek — Cc). See Appendix C for the actual expression of 

We have determined the critical Coulomb interaction Ud solving eq. (46). For a hypercubic 
lattice, we obtained [/ci(MLA) = 3.237, while we find C/ci(HF) = 3.693 when we adopt the 
Hartree-Fock wave function. The reduction of Ud due to electron correlations is understood 
from the DOS in Fig. 6. There we observe that the electron correlations on the static matrix 
elements enhance the Mott-Hubbard peaks at a; = ±U/2. It reduces the amplitude of DOS 
near the Fermi level, and therefore lowers the critical value Ud- 

The present value [/ci(MLA) = 3.237 quantitatively agrees with the numerical result 
[/ci(NRG) = 3.25 obtained by the NRG.^o In the case of the Bethe lattice for which the 
noninteracting DOS is given by p{e) = -k'^^JI — e^, we obtained ?7ci(MLA) = 3.359 and 
[/ci(HF) = 3.827. The critical value for the Bethe lattice has been calculated by various 
methods. These results are summarized in Table I together with the present result 
(MLA+RPT-0) . There are some discrepancies in Ud among the theories in case of the Bethe 
lattice. Our result agrees well with Ud = 3.39 obtained by the projective self-consistent 
technique (PSCT)^^ which is exact in the low energy region, and Ud = 3.41 obtained by 
the local moment approach (LMA).^^ Note that the results of the ED (Exact Diagonalization 
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Method)^^ and NRG^'^ methods considerably depend on the way of Hne-broadening for the 5- 
function spectrum in the case of the Bethe lattice. The present approach describes Ud within 
1% error for both the hypercubic and Bethe lattices. 

The quantitative description of the excitation spectra for the Mott insulator with use of 
the RPT-0 may come as a surprise because the RPT is an approach starting from the weak 
Coulomb interaction limit due to the expansion of the screened memory function with respect 
to the interaction Liouvillean.*^ In this respect, it is worth pointing out that the RPT-0 can 
also describe the dynamics of the strongly correlated electrons via the effective medium and 
atomic self-energy in the denominator in eq. (41) since the RPT-0 memory function reduces 
to the Hubbard III approximation^^ in this region. 

5. Summary 

We have proposed a new local-ansatz wavefunction with momentum-dependent variational 
parameters (ML A) to improve the LA by StollhofF and Fulde.^ It is constructed by using the 
'flexible' local operators which produce the two-particle excited states in the momentum space 
from the Hartree-Fock state and project those states onto the local excited states in the real 
space. The best wavefunction is chosen by controlling the momentum dependent variational 
parameters of the excited states in the momentum space on the basis of the variational prin- 
ciple. We obtained the ground-state energy of the MLA within a single-site approximation. 
Minimizing the energy, we derived a self-consistent equation for the variational parameters, 
and obtained an approximate solution which interpolates between the weak Coulomb interac- 
tion limit and the atomic limit. The correlation energy in the MLA agrees with the result of 
the second-order perturbation theory in infinite dimensions in the weak Coulomb interaction 
limit and yields the correct atomic limit as it should be. 

We have investigated numerically the validity of the theory using the half-filled band 
Hubbard model in infinite dimensions. We verified that the MLA improves the LA in the 
whole range of the Coulomb interaction energy parameter U. For the hypercubic lattice, we 
found that the MLA yields the correlation energy lower than that of the GA in the range 
< U < 3.28. The double occupation number in the MLA is smaller than that of the LA 
irrespective of U. The GA overestimates the double occupancy in the range < ?7 < 2, and 
underestimate it in the range U > 3. We found that the MLA shows a reasonable energy 
dependence of the momentum distribution in the range < U < 3.0. This is qualitatively 
different from the LA and the GA because both of them lead to the energy-independent 
momentum distributions below and above the Fermi level. The quasiparticle weights in the 
MLA are very close to those of the NRG in the range < U < 2.5, while the LA and the 
GA overestimate them in general. These results suggest that the MLA is applicable to the 
systems with U/W < 1.5, for example, the systems like transition metals and alloys. Here W 
denotes the band width of the noninteracting system. 
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The critical Coulomb interaction Uc2 in the MLA was obtained from the vanishment of 
the quasiparticle weights as Uc2 = 3.21. It is comparable to C/c2 = 4.10 in the NRG, while the 
LA and the GA give larger values Uc2 = 7.82 and 4.51, respectively. 

We have also shown that the MLA combined with the PM-CPA is useful for understanding 
the correlation effects on the excitation spectra in the insulator regime. The MLA wavefunction 
allows us to calculate the static correlations in the retarded Green function obtained by the 
PM-CPA. We have demonstrated that the MLA+PM-CPA can quantitatively describe the 
excitation spectra in the insulator regime. Calculated critical Coulomb interactions Ud for a 
gap formation agree with the best results obtained by the other methods within 1 % error for 
both the hypercubic and the Bethe lattices. 

Although the present approach interpolates between the weak Coulomb interaction limit 
and the atomic limit, and much improves the LA, it does not describe the metal-insulator 
transition in a self-consistent way. The wavefunction continuously changes from the Hartree- 
Fock metallic state to the atomic one in the present theory, and does not show any anomaly 
at Uc2 {Uci) obtained from the momentum distribution (the excitation spectra with use of the 
MLA+PM-CPA). Further improvements of the theory toward the strongly correlated region 
are desired to describe the metal-insulator transition in a self-consistent way. 
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Appendix A: Average {A ) in the single-site approximation 

In this Appendix, we derive the formula (13) in the single-site approximation. Let us 
consider Aj^ and Sjv such that 



An = - d\)\A[\{{l - O.. 



(A-l) 



BN = {[\[{l-Ol)\[\[{l-d,) 



Expanding Bjy with respect to site 1, we obtain 

(1)., ^txir-n-a), 



B 



N 



[n'"( 



i-oj 



Oi 



l-Oi 



and 



«2^=([n'"(i-oi)][n'"(i-« 



(1), 



(A-2) 



, (A.3) 



(A-4) 
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Here the product II/^^ means the product with respect to all the sites except site 1. 

When we calculate Bj^ applying Wick's theorem, we neglect the contractions between 
different sites. This is a single-site approximation and then eq. (A-3) is expressed as 

Bn = {{i-d\){i-d,))yZi- (A-5) 

We can make the same calculations for Ajq. In this case, there are two-types of terms, the 
terms in which the operator 0\ is contracted to A and the other terms with 0\ contracted to 
the operators Oi {i^l). We have then 

An = <(l-Ol)i(l-Oi)>oi?Ei + <(l-Ol)(l-Oi)>34)_, , (A.6) 

and 

4'-i=([n"'(i-oi)]-i[n"V-o<)])„- 

i i 

Successive application of the recursive relations (A-5) and (A-6) leads to 

An = E<(l-Ol)^(l-00>o^S^-i' (A-8) 

i 

Bn = {{1- Ot) (1 _ 0.)>, B^l, = n ((1 - Oj) (1 - 0,))o ■ (A-9) 

i 

Taking the ratio An/Bn, we reach eq. (13). 

(A) = y ^—S . (A-IO) 

Appendix B: Laplace transform for the correlation calculations 

Using the relation 

POO 

= -i die*(^-"4+e3-e2+ei+.c)t ^ 

Jo 



(B-l) 



z — e4 -I- €3 — 62 + ei + ec 
we can reduce the number of integrals in the physical quantities. Here z = lo + id, and S 
is an infinitesimal positive number. Laplace transforms of various elements in the physical 
quantities are summarized as follows. 

{Hdi)o = {d\H)i 

POD 

= if/2[l-,7(l-2(n,t)o)(l-2(ni;)o)] / dte''^' a^{-t)ai{-t)b^{t)bi{t) . {B-2) 

Jo 

{0\Hdi)o = -C/2[l-??(l-2(n,t)o)(l-2(n^)o)]2 

/•CO 

X / dtdt'e''^^^+^'\a^{-t-t')h-^{t + t')ai{-t-t')hi^{t + t') 
Jo 

-at(-i - t%{t + t')ai^{-t - f')bi{t + t') 
+a^i-t - t')b^{t + t')aii-t - t')bi{t + t') 
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-ait(-t - t%{t + t')aii-t - t')hi{t + t')] , (B-3) 

{dlOidi)o = -[/2[l-r?(l-2(n,t)o)(l-2(n,i)o)P 

/•oo 

X / dtdt'e''^^^+^'^[a^{-t)b^{t + t')ai{-t)bi{t + t')a^{-t')ai{-t') 
Jo 

-a^i-t)b^it + t')ai{-t - t')biit)a^i-t')bi{t') 
-a^i-t - t%it)aii-t)b^{t + t%it')a^i-t') 
+a^^-t-t%{t)a^{-t-t%{t)b^{t')b^{t')] , (B.4) 

(O]"a)o = -i72[l-r/(l-2(n,T)o)(l-2(n,i)o)]' 

POD 

X dtdt'e''^^^+^''^a^{-t-t%{t + t')ai{-t-t')bi{t + t') . (B-5) 



Here 



aa{t) = J dep{e)f{e + e<,) e"^^* , (B-6) 
ba{t) = J dep(e)(l - /(e + e,)) e"*^* , (B-7) 
aia{t) = [ dep{e)f{e + e^) ee"^^* , (B-8) 



X 







X 



bia{t) = J dep{e){l-f{e + ea))ee-''' , (B-9) 

and Co- = eo + U{ni-a-)o — M is the Hartree-Fock level measured from the Fermi level p. p{e) 
in the above expressions denotes the density of states for e^, i.e., the Fourier transform of tij. 
Correlation contribution to the momentum distribution function (36) is given by 

N{dlhkadi)o = -2C/2[l-r/(l-2(niT)o)(l-2(n,|)o)]' 

dtdt'e'^^^^+^'^aait + t')bai-t - t') 

[f{eka)a-a{t + t')e-^^^(*+*') - (1 - /(efc,))a_,(-t - t')e'^^ ■ (B-10) 

Correlation contribution to the electron number (35) which appears in the calculation of the 
double occupation number is expressed as 

{0\hM)Q = -C/2[l-r7(l-2(nit)o)(l-2(n^)o)]2 

X / dtdt'e'''=^*+^'^[a-a{-t-t')b-„{t + t')a„{-t-t')b^{t)ba{t') 
Jo 

-a-a{-t-t')b-a{t + t')aa{-t)b^{t + t')aa{t')] . (B-11) 
The element (37) for the calculation of the double occupancy is expressed as 

{OlOi)o + {OA)o = -2iU[l - 7j{l - 2(ni|)o)(l - 2{ni^)o)] 



20/23 



J. Phys. Soc. Jpn. Full Paper 

POO 

X / dte'''^'a^i-t)b^{t)a^{-t)biit) . (B-12) 
Jo 



Appendix C: Expressions of X{e,e',e"), GqI{z), and c^"^ in the MLA 

In the R = approximation, the correlation function X[e,e',e") given by eq. (44) is 
obtained from the formula (13) as 

X,ie,e',e") ^ X2{e,e',e") 

1 + (OjOi)o 1 + {dldi)o 

For the half-filled band, we have 



Xie, e', e") = x{e, e', e") - ^^i^^^^ + > . (C-l) 



X(e, 6', e") = /(-e)/(-eO/(e") + /(6)/(e')/(-e") , (C-2) 

Xi(e,e',e") = -f{^')f{-e)v{e-^')-f{-^')f{e)v{^'-e) 

+f{-e')fi-eMe + e') + f{e')f{e)u{-e - e') , (C-3) 

X2ie,e',e") = f{-e')f{e")K^{e'-e")-f{e')f{e")K2{-e\e") 

-f{-e')f{-e")K2{e', -e") + f{e')f{-e")K,{e" - e') 
-f{e")f{e)n2i-e", e) + /(e")/(-e)Ki(e - e") 
+/(-e")/(e)Ki(e" - e) - f{~e")f{-e)K2{e\ -e) 
+f{e')f{e)Ki{-e' - e) - f{e')f{-e)K2{e,e') 
-f{-e')f{e)K2{e', e) + /(-e')/(-e)Ki(e' + e) 
+{f(e)-f{e"mf{-e')\{e')-fie')X{-e')] 
-ifi-e) - f{e'))[f{e")\{-e") - fi-e")X{e")] 

+{f{e") - f{e'))[f{e)X{-e) - f{-e)X{e)] . (0-4) 

Here u{e) is defined by 

z/(efe) = J dede'p{e)p{e')f{e)f{-e'W - e + e^) , (C-5) 

•q{e) is defined by r7(e) = U/{e — ec). /«i(e), K2(e, e'), and A(e) in X2 are defined by eqs. (49), 
(50), and (51), respectively. 

Substituting eq. (C-l) into eq. (43) and making use of the Laplace transform (B-1), we 
obtain the expression for gS(z) with use of the Laplace transform as 

G£(.)=M„(.-E(.)) + iH2(__|£)). (C.0) 

1 -I- {U-Ui)o 

Mo{z) = -i j dte'^\b{-tf + b{tf) , (C-7) 
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M2{z) = -3i J dte''^(P{t)[Ki{t) - K2it) + {b{-t) - b{t))Xi{t)] , (C-8) 

ni{t) = -U^ J dt'dt"e''^^''+'"\b{t + t'f + b{t' - tf)b{t' + t"fb{t"f , (C-9) 

. -2U^ J *'*"e-<"-).(t>(t + mfMf - tW + f? , (CIO) 

Ai(t) = j dt'dt"e'''^^*'+^"\b{t + 1') - b{t' - t))b{t' + t"fb{t") . (C-11) 

Here ^(i) = J dep(e)e'^*, and we used the relation a(t) = 6(— t) for the half-filled and symmetric 
band. 

(2) 

In the same way, the correlation contribution C2 to the second moment of the memory 
function (i.e., eq. (48)) is expressed as 

= -12U^ j dtdt'e''^^'+''h{t + t'fbi{t)b{t') [bi{t)b{t') + b{t)bi{t') -ab{t + t')] . (C-12) 

Here a = l/\/7r for the hypercubic lattice and a = 4\/2/37r for the Bethe lattice. b{t) and 
bi{t) are defined by eqs. (B-7) and (B-9), respectively. 
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